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Abstract 

o 

■ The point process of vertices of an iteration infinitely divisible or more specifically 

^ . of an iteration stable random tessellation in the Euclidean plane is considered. We 

explicitly determine its covariance measure and its pair-correlation function as well 
as the cross-covariance measure and the cross-correlation function of the vertex point 
process and the random length measure in the general non-stationary regime, and 
we give special formulas in the stationary and isotropic setting. Exact formulas 
^ ■ are given for vertex count variances in compact and convex sampling windows and 

asymptotic relations are derived. Our results are then compared with those for a 
Poisson line tessellation having the same length density parameter. Moreover, a 
functional central limit theorem for the joint process of suitably rescaled total edge 
count and edge length is established with the process t > 0, arising in the 

limit, where £ is a centered Gaussian variable with explicitly known variance. 
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1 Introduction 

Random tessellations have attracted particular interest in stochastic geometry because of 
their wide applications ranging from classical geological problems to recent developments 
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in telecommunication, see [10] and [13]. It is one of the main purposes of the related theory 
to develop new classes of random tessellations that are mathematically tractable and yet 
rich enough in structure so that they may serve as new reference models for applications 
beside the classical Poisson hyperplane and the Poisson-Voronoi tessellation. A very recent 
model, the so-called random STIT tessellations - the abbreviation STIT stands for stable 
under iteration, see below - was introduced by W. Nagel and V. WeiB in [8]. One of 
the main features of these tessellations that distinguishes them from the above mentioned 
model classes is the property that their cells are not side-to-side, see Figure 1. This causes 
new geometric effects whose planar first-order properties in terms of mean values were 
explored in [9]. 

This paper will deal with second-order characteristics and central limit theory for the pla- 
nar case, a topic that has been considered first by WeiB, Ohser and Nagel in [16] and 
later also studied in [12], where beside other characteristics the variance of the total edge 
length was determined. Here, in contrast, we will deal with second-order properties of 
the point process of vertices of the tessellation. To provide more general results, we will 
not restrict here to the class of stationary random iteration stable tessellations (STIT tes- 
sellations), and instead we will study the larger class of non-stationary random iteration 
infinitely divisible tessellations in the plane. They were introduced in [12] as generaliza- 
tions of STIT tessellations and in Section 2.1 we will recall their construction in the spirit 
of Mecke-Nagel- WeiB (MNW) [8] . It is an important observation that this spatio-temporal 
MNW-construction can be interpreted as a continuous time Markov process on the space of 
tessellations, whence the general theory of Markov processes is available. Before extending 
some mean value relations from the stationary iteration stable to the non- stationary itera- 
tion infinitely divisible case in Section 2.3, we will formulate the main technical tools from 
the theory of Markov processes in Section 2.2, on which our main results are based. They 
are the content of Sections 3 and 4. The variance of the total number of vertices and that 
of the total number of maximal edges in a bounded observation window will be calculated 
in Section 3.1 for very general driving measures, whereas in Section 3.4 we specialize to 
the motion-invariant case. The vertex pair-correlation measure for general line measures 
is considered in Section 3.2 and the considerations are specialized again to the stationary 
and isotropic setting in Section 3.5. Moreover, we determine in Section 3.3 the exact cross- 
covariance measure of the vertex point process and the length measure in the general case 
and provide specialized formulas in Section 3.6 for the stationary and isotropic regime. 
Here, the most explicit formulas are available. Another topic treated there concerns the 
variance asymptotics for a sequence Wr of growing windows as R tends to infinity. Based 
on these results, Section 4 deals with the central limit problem. It will be shown that a 
certain rescaled bivariate process of edge counts and edge lengths of a stationary random 
iteration stable (STIT) tessellation converges to the process (£,££)) where £ is a centered 
normal random variable with an explicitly known variance. 

We would like to emphasize that our results may be of interest in the context of statisti- 
cal model fitting of random tessellations to real data, see [10]. Our central limit theory 
could be a base for statistical inference of tessellation models and related functionals, and 
asymptotic confidence intervals and statistical tests can be derived from them, since we 
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Figure 1: Realizations two STIT tessellations with different driving measures, the invari- 
ant measure (left) and a measure concentrated on lines pointing in only two orthogonal 
directions with weight 1/2 (right) 



make the first- and second-order moments explicitly available. 

For general notions and notation related to stochastic geometry and the theory of random 
tessellations, which are used in this paper, we refer the reader to [11] and [13]. 

2 Background Material 

2.1 Iteration Infinitely Divisible and Iteration Stable Random 
Tessellations 

Consider a compact and convex set W C M 2 with non-empty interior and a diffuse (non- 
atomic) measure A on the space of [M 2 ] of lines in the plane enjoying the local finiteness 
property, stating that A([C]) < oo for any compact C C M 2 , under the usual notation 

[C] = {lines L : L n C ^ 0}, 

letting [C] stand for the set of lines that have non-empty intersection with C. Below we 
briefly describe the construction of an iteration infinitely divisible tessellation in W, the 
construction is called the MNW-construction to honor its inventors - Mecke, Nagel and 
WeiB , who introduced it in [8]. To begin, assign to W an exponentially distributed random 
lifetime with parameter given by A([W^]). Upon expiry of this random time the cell W dies 
and splits into two sub-cells W + and W~ separated by a line in [W] chosen according to the 
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law A([VF]) _1 A(- fl [W]). The resulting new cells W + and W~ are again assigned indepen- 
dent exponential lifetimes with respective parameters A([W + ]) and AQV4 7- ]), whereupon 
the construction continues recursively and independently in each of the sub-cells W + and 
W~ until some deterministic time threshold t > is reached. The random tessellation con- 
structed in W until the time t will be denoted by Y(t, W). To put it formally by Y(t, W) 
we understand here the random closed subset of W arising as the union of the boundaries 
of cells constructed by the time t. The cell-splitting edges are called the maximal edges 
(in the related literature often called I-segments, as assuming shapes similar to that of 
the literal I) of the tessellation Y(t, W) and the family of all such edges is denoted by 
MaxEdges(F(t, W)). Note that such maximal edges can get further subdivided between 
their birth time and the time t, that is to say in the course of the MNW dynamics there can 
appear additional vertices in their relative interiors. Thus, we have a distinction between 
maximal edges with a possible interior structure and those edges which are not maximal 
(for example the sides of cells or the primitive elements that are bounded by vertices but 
have no interior structure). 

In general, the random tessellations Y(t, W) in W do not have to arise as windowed re- 
strictions of stationary (stochastically translation invariant) or isotropic (stochastically 
rotation invariant) processes. If we assume in addition though that the driving measure 
A has, under polar parametrization, the product structure A = r£ + (g) 1Z with r a positive 
constant, £ + the Lebesgue measure on M + and a spherical directional distribution 1Z on 
the unit circle «Si, then Y(t, W) has the property of being a suitable windowed restriction 
of a stationary random field. The random tessellation Y(t, W) is furthermore isotropic in 
a similar restriction-wise sense if and only if 1Z is the uniform distribution, i.e. if and only 
if A is a multiple of the standard motion-invariant measure A iso := £ + <g> v x of unit length 
intensity on the space of lines in the plane, where v\ is the uniform distribution on S\. 
We review now some of the important properties of the tessellations Y(t,W), proofs of 
which can be found in [7], [8] and [12]. We formulate them only for the planar case even 
if they are true in higher dimensions as well. 

1. The random tessellations Y(t, W) are consistent in W, that is to say Y(t, V) fl W = 

Y(t, W) for W C V for V, W compact convex, where = means equality in distribution. 

This implies the existence of the whole- plane tessellation Y(t) such that Y(t, W) = 
Y(t) fl W for each compact convex W. 

2. The MNW-construction satisfies the Markov property with respect to iteration of 
tessellations in the time parameter t, i.e 

Y(s, W) ffl Y(t, W) = Y(s + t, W), 

where ffl denotes the operation of iteration for whose exact definition we refer to 
the above mentioned papers. This operation can roughly be explained as follows: 
Let Yq be a random tessellation and let Ti, T2, . . . be a sequence of independent and 
identically distributed random tessellations in the plane. To each cell of the frame 
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tessellation Y we associate the tessellation T^. Now, we make a local superposition 
of Yq and the cutouts of Tjt in each cell of the frame tessellation. If we let Y 1 be 
a tessellation with the same distribution as Tk for any k = 1,2, ... , we denote the 
outcome of this procedure of local superposition by Y 53 Y\ and call it the iteration 
of Y with Fl 

3. It directly follows from the above Markov property that the random tessellations 
Y(t, W) are infinitely divisible with respect to iteration, i.e. 

Y(t, W) = Yit/n, W) EH ... EE Yit/n, W) 

V v ' 

n 

for any nGN and any compact convex W . For this reason, the tessellations Y(t, W) 
are referred to as iteration infinitely divisible random tessellations. However, it is 
worth pointing out that it is currently not known whether any iteration infinitely 
divisible random tessellation can be obtained by the MNW-construction. If in addi- 
tion the driving measure A is translation invariant, it can be verified that Y(t, W) is 
even stable under iteration, i.e. 

Y(t) = n( Y(t) EH . . . mY(t)) , n G N, 

n 

where n(-) means the dilation with a factor n, i.e. nY = {n ■ x : x G Y}. Note 
that this equation must be understood symbolically, since properly the operation 
of iteration involves a sequence of independent and identically distributed random 
tessellations, but we adopt here the usual convention from the earlier work [7]. Thus, 
in the stationary case we are in fact dealing with iteration stable random tessellations 
or STIT-tessellations for short. 

4. The intersection of an iteration infinitely divisible random tessellation Y(t, W) having 
driving measure A with an arbitrary line L G [W] is a Poisson point process with 
intensity measure A h-> tA(L4]) with Ac W fl L a Borel set. In particular, for 
x,y eW the probability that they belong to the same cell of Y(t, W) equals 

F(x and y are in the same cell of Y(t, W)) = e ~ tA{ ^\ 

where by xy C W we mean the line segment joining the points x and y. 

2.2 Martingales in the MNW Construction 

As already seen in the introduction above, the MNW-construction of iteration infinitely 
divisible or iteration stable random tessellations Y(t, W) in finite volumes W C M 2 for 
general locally finite and diffuse driving measures A enjoys the Markov property in the 
continuous time parameter t. In our previous work [12] we have used this fact combined with 
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the classical theory of martingale problems for pure jump Markov processes to construct 
a class of natural martingales associated to the MNW process. In this paper we only need 
a part of that theory. To formulate the required results, for a bounded and measurable 
functional of a line segment (tessellation edge) and for a tessellation Y, usually taken to 
be a realization of Y(t, W) for some t > 0, write 

MY) = E ^ e )- 

eGMaxEdges(y) 

Note, that with Y as above, for each straight line L G [W] the intersection Y H L is 
just a one-dimensional tessellation of L HW which can be identified with the collection 
Segments (Y n L) of its constituent segments. Bearing this in mind we write 

MY) = f E <Ke)k{dL). (1) 

It is also convenient to introduce the bar notation for centered versions of these quantities 
with Y = Y(t, W), that is to say 

E*(y(t, W)) := ^(Y(t, W)) - EE (y(t, W^)) 

and likewise 

A^(Y(t, W)) := A^(Y(t, W)) - EA^Yit, W)). 
With this notation, in view of (41) in [12] we have 

Proposition 1 For bounded measurable segment junctionals and ip the stochastic pro- 
cesses 

^(Y(t,W))- f A^Y(s,W))ds (2) 
Jo 

and 

E4(Y(t,W))Ei,(Y(t,W)) - [ A H (Y(s,W))ds- 

Jo 

- [ [At(Y(s,W))^(Y(s,W))+A^(Y(s,W))t:^Y(s,W))}ds (3) 
Jo 

are martingales with respect to the filtration $$ t generated by (Y(s, W))o< s < t . 

2.3 Mean Values for Edge and Vertex Processes 

This subsection recapitulates some basic first order properties of iteration infinitely divisible 
tessellations, mostly known in the stationary set-up, for sake of reference in further sections. 
Using (2) for <p = 1 yields, upon taking expectations, 

EEi(y(t, W)) = [ EA 1 (Y(s, W))ds. (4) 
Jo 
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However, in view of the definition (1) it is easily verified that 

A^Y) = A([W]) + A([ e ])=A([W]) + E AH (y), (5) 

eeMaxEdges(Y) 

where A[-] stands for the edge functional e h-> A([e]). Applying (2) once again for = A[-] 
leads to ^ 

EZ A[ . ] (Y(t,W))= [ EA A[ .](y(s,W))ds. 
./o 

But 

A A[ . ] (y(s,W))= / A([LnW])A(dL) = ((AnA))(W), Vse(0,t], (6) 

where the locally finite point-intersection measure ((ADA)) on IR 2 is given by 

((AnA»:=/ / 6 LnL ,A(dL>),A(dL) (7) 

J[R 2 ] J[L] 

where, recall, 5t.\ stands for the Dirac unit mass at the argument, so that in other words 
((A n A)) (A) = (A x A){(Lx, L 2 ) e [A] x [A], L x n L 2 e A}, A C R 2 . 

Hence 

ES A[ .](y(s,W))=s((AnA»(W). (8) 
Combining (4), (5) and (8) finally yields 

EE^y^, W)) = tA([H/]) + ^((An A»(W). (9) 

Note, that the formula (9) when specialized to the translation-invariant set-up contains an 
extra boundary correction term iA([W]) in comparison to the classical mean value formula 
for the iteration stable (STIT) random tessellations as given in [9], which says that the 
density of maximal edges in W is just ^-((AnA))(W). This additional boundary correction 
term iA([W]) comes from the fact that we count edges rather than edge midpoints. Thus, 
it can happen that in two neighboring regions one observes two distinct edges which may 
coalesce into one edge when putting these regions together into one area. 

3 Second Order Theory for Edge and Vertex Pro- 
cesses 

3.1 Variance Calculation for the General Case 

We consider the most general case first and study iteration infinitely divisible random 
tessellations Y(t, W) with general locally finite and non-atomic driving measures A. We fix 
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t > and a compact and convex observation window W C M 2 as in Section 2.1. 
First we use (3) with <j> — tp = 1 to conclude that 

El(Y(t,W))- f A 1 (Y(s,W))ds-2 f A l (Y(s,W))i: 1 (Y(s,W))ds 
Jo Jo 

is a martingale with respect to 9^. Using (5) and taking expectations we get 

Vai^Y^W))) = tA([W})+ [' EE m (Y(s,W))ds 

Jo 

+2 / Cov(EA H (y(s,WO),Ei(y(s,WO))<fe. 
Jo 

It remains to find an expression for the covariance Cov(£ A [.](F(s, W)), T, 1 (Y(s : W))). Such 
an expression can be found by applying (3) once again, this time with = 1 and i[> — A[-] 
and t replaced by s and s by u. We obtain in this way 

Cov(E AN (y( s ,wO),£i(ns,wO)) = / EA A[ .](y(«,wO)du 

Jo 

+ f Cov(A 1 (y(^, iy)), E A[ .](y(u, wO))d« 

./o 



(10) 



+ / Cov(A A[ . ] (y(u,w)),i; 1 (y(u,w)))d«. 

In view of (6), A A u(-) is a constant and hence the covariance involving it vanishes. Resort- 
ing again to (5) we end up with 

Cov(£ A [.](y(s, W)), E!(y(s, W))) 

= [ EA A[ . ] (Y(u,W))du+ [ Var(E AN (y(u,W)))^. (11) 
Jo Jo 

Putting together (10) with (11) yields the following expression for Var(Ei(y(£, W))): 

Vax(E 1 (Y(t,W)))=tA(\W])+ [ EE A [.](y(s, W))ds 

Jo 

+ 2[ f [ EA A[ . ] (Y(u,W))duds+ f [ Var(E AN (y(u, W))duds ) . (12) 
\Jo Jo Jo Jo J 

It remains to find Var(£ A [.](y(t, W))). To find an expression, apply (3) again with <fi = 
ip = A[-] to get, upon taking expectations and using that A A [.] is a constant as remarked 
in (6) above, 

Var(£ A[ .](y(t,W0))= / EA {A[W (Y(s,W))ds. (13) 

Jo 
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However, using (1) 

A (A [.])a(y(s, W)) = / / / 1[L n Li and L n L 2 are 
V[w] ./[.Lnw] ./[Lnw] 

in the same cell of W)]K(dL 1 )K(dL 2 )K(dL) 

and hence 

EA (A[ . ])2 (y( s ,w)) 

/ / P(L n Li and L f) L 2 are 
[W] ./[Lnw] ./[Lnw] 

in the same cell of F(s, W0)A(dLi)A(dL 2 )A(dL) 



/ 



/ / exp(-sA([L(L 1 ,L 2 )]))A(dL 1 )A(dL 2 )A(dL) (14) 



[W] J [Lnw] J [Lnw] 

where L(Li,L 2 ) stands for the segment joining the points L D L\ and L fl L 2 and where 
the last equality follows by the property 4 of the tessellation F(t, W) as listed in Section 
2.1. 

To neatly formulate our theory, denote by (((A x A) fl A)) the segment-intersection measure 
on the space [M 2 ] of finite linear segments in R 2 given by 

(((A x A) n A)) = / / / 5 L(Ll , L2) A(rfL 1 )A(rfL 2 )A(rfL) (15) 

J[R 2 ] J[L] J[L] 

and observe that this defines a locally finite measure, charging finite mass on collections 
[A] of segments with both ends falling into a bounded set A C R 2 , because of the local 
finiteness of A. With this notation, combining (6), (8), (11), (12), (13) and (14) yields the 
following result: 

Theorem 1 For general locally finite and diffuse driving measures A, denoting by 



n-l 



k=n fc=0 

the n-th tail of the exponential series at u, we have 

Tr(-tA([e])) 



Var(E AH (y(f,W))) = - / Tl ^ (N)) ({(AxA)nA))(rfe) 

J[w] A ([eJ) 



and 



Cov(E A[ . ] (y(t,^)),i: 1 (y(t,w))) 

^ xp (-tA([e]) ; 

'[w] 



= t((A n A»(wo + J- w Tr \ll { } e])) i(( A x A ) n A »( rfe ) 
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and 

3f 2 

Var^Y^WO)) = tA([W\) + ^«A n A))(W) 

- 2 l ^T^^ (((AxA)nA))(de). (16) 
J[w] A (Nr 

3.2 Vertex Pair-Correlations for the General Case 

Also in this subsection we stay in the general set-up of locally finite and diffuse A. We shall 
extend here the calculations made in Subsection 3.1 above to determine the pair-correlation 
structure of the vertex point process Vy(t,w) generated by Y(t, W). For definiteness we 
adopt the convention that Vy(t,w) does not include the boundary vertices, this way each 
vertex arises at intersection of exactly two maximal edges. Recalling the consistency rela- 
tion Y (t, W) = Y(t) n W, we see that the covariance structure between bounded regions 
U, V C M 2 does not depend on W as soon as both U and V are contained in the interior of 
W. To put this in formal terms, consider the whole-plane covariance measure Cov(Vy(t)) of 
the point process Vy(t) on (IR 2 ) 2 = K 2 xt 2 (also called the second-order cumulant measure) 
given by the relation 

/ (/® g)dCov(V n t)) = Cov(Z vf (Y(t)),Z vg (Y(t))) (17) 

J(R 2 ) 2 

holding for all /, g : 1R 2 — > K. bounded measurable and of bounded support, where r/ is the 
edge functional 

V f (e) = J2 /(*) 

:r€Vertices(e) 

and likewise for rf . Note, that, even though we are apparently dealing with functionals 
T, v f(Y(t)) and S^ 9 (F(t)) defined on the whole-plane process, they can be safely replaced 
by T* v f(Y(t, W)) and T* v g(Y(t, W)), respectively, for some W containing the supports of 
/ and g, hence our martingale relations given in Proposition 1 hold here with no extra 
assumptions. For the same reasons all integrals below with apparently unbounded inte- 
gration domains are effectively bounded due to the bounded supports of / and g, which 
we are going to exploit without further mention. It is readily seen from (1) that for each 
possible realization Y of Y(t) or Y(t, W) in a domain W containing the supports of /, 

A ¥ (Y) = 2 f J2 f( x ) A (dL) = 2£ A/[ . ] (y), (18) 

where A^[e] = f(e fl L)A(dL) and the factor 2 comes from the fact that each point of 
the tessellation is contained in exactly two maximal edges. Consequently, using (3) for 
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(j) = rjf, ip = rf and taking expectations, we get 

Cov(E T? / (Y(t)), E^ g (y(i))) = f\A vSr , g (Y(s))ds 

Jo 

+2 / Cov(E A/[ . ] (y( s )),E, 9 (y( s )))d s 

Jo 

+2 / Cov(E As[ . ] (y( S )),E r , / (y( S )))d S . (19) 

Proceeding as in the previous Subsection 3.1, we turn now to the calculation of the covari- 
ance Cov(E A /[.](Y(s)), E T?9 (y(s))). To this end we note that 

A A f [ . ] (Y)= [ Af[L]A(dL)= [ M(AnA)), (20) 

J[«. 2 ] Jr 2 

whence A A f^(-) is a constant, and we use again (3) for <fi = A* , tp — rj 9 to get in view of 
(18) 

Cov(E A/[ . ] (F( S )),E, 9 (r( S ))) 

= [ EA Af[ . ]vB (Y(u))du + 2 [ Cov(E A / [ . ] (y(u)),E A9 [.](y(«)))d«. (21) 
Jo Jo 

Finally, one further use of (3) with = , ip = A 9 and application of (20) yields 

Cov(E A/[ . ] (y( M )),s A9[ . ] (yH))= / EA A/[ . ]A9N (y(t;))^. (22) 

Jo 

Thus, using (19), then twice (21), once with / and g interchanged, and then (22), we get 



Cov(E„/(y(*)),£^(y(t))) = / EA ¥va (Y(s))ds 



o 



+2 EA Af[ . ]vg+Aa[ . ]rjf (Y(u))du + 8 / / / EA^^^Y (v))dvduds, 

Jo Jo Jo Jo Jo ^23) 

since A A f^ v9 + A Ag ^ v f = ^4 A /[.]^s+ A s[.]r,/- It remains to calculate the expectations of the 
functionals present in these integrals. However, this is easily done by recalling that, 
for L G [1R 2 ], the intersection Y(t) H L is the Poisson point process with intensity measure 
L D A !->■ £A(L4]), see property 4 in Subsection 2.1, whence, in view of (1), 



EA+(Y(t)) = \l 0(e)exp(-tA([e]))(((tAxtA)nA))(rfe) 

'[9?] 



0(e) exp(-tA([e]))(((A x A) n A))(de) (24) 
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for bounded measurable such that 0(0) = and locally defined in the sense that there 
exists a bounded convex W such that 0(e) = 0(e D W) for all e. Note that the extra 
prefactor of 1/2 comes from the fact that the segment-intersection measure defined by (15) 
counts each segment twice, once for each of the two orderings of its two termini. Putting 
(24) together with (23) yields now 

Cov(S r , / (F(t)),S r , 9 (F(t))) = 

+ \ \ V f(e)r 1 %e)l 1 (s 2 eM-tH[em) (((Ax A)nA»(de) 

1 J[R 2 ] 

+ I _(A f [e]rf[e] + A 9 [e]rj f (e)) X 2 {s 2 exp(-sA([e])); t) (((A x A) nA))(de) 



+ 4/ A / [e]A fl [e] X 3 (s 2 exp(-sA([e])); t) (((A x A) n A))(de), (25) 

J[R 2 ] 

where the multiple integral X n is given by 

ft rsi rsn-i 1 rt 

l n (f(s);t):= / .../ f(s)dsds n _ 1 ...ds 1 = - — (t - s)"" 1 f(s)ds 

Jo Jo Jo \ n ~ L)\ j Q 

so that in particular 

X 1 (s 2 exp(-As);t) = A~ 3 (2 - (A 2 t 2 + 2Xt + 2) exp(-At)), (26) 

X 2 (s 2 exp(-As);t) = A" 4 (2At - 6 + (A 2 t 2 + 4At + 6) exp(-At)), (27) 

X 3 (s 2 exp(-As);t) = A~ 5 (A 2 t 2 - 6 At + 12 - (A 2 t 2 + 6\t + 12) exp(-Af)). (28) 

For a segment (edge) e consider the measures A e and A[- n e] on R 2 that are defined by 

A e := £ S * 

xeVcrtices(e) 

and 

(A[-ne])(A) = A([Ane]), A C R 2 . 

With this notation, putting together (17) and (25) yields in view of the definitions of r\' 
and A-* [•] the following 

Theorem 2 For general locally finite and diffuse driving measures A we have 

Cov(Vy W )= / ^(A e ®A e )X 1 ( S 2 exp(- S A([e]));t)(((Ax A)nA))(de) + 

J[R 2 ] Z 

I (A e <g> A[- n e] + A[- He]® A e ) X 2 (s 2 exp(-sA([e])); t)(((A xA)n A»(de) + 

4/ (A[-ne]®A[-ne])X 3 (s 2 exp(-sA([e]));t)(((Ax A)nA))(de). (29) 
JTR 2 ! 
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An intuitive understanding of the structure of the covariance measure in Theorem 2 comes 
by noting that the first integral in (29) takes into account pairs of vertices constituting 
ends of the same maximal edge, the second one corresponds to pairs of vertices with the 
property that one of them is an internal vertex of a maximal edge of which the second 
point is a terminus, whereas the third integral corresponds to pairs of vertices lying on the 
same maximal edge but not being its termini. Thus, other pairs of vertices (not lying on 
the same maximal edge) are not present in the covariance measure, roughly speaking this 
is because the maximal edges are the only means of propagating dependencies in iteration 
infinitely divisible tessellations, an intuition to be made more specific in our further work 
in progress. 

3.3 Edge- Vertex Correlations in the General Case 

In this subsection, still placing ourselves in the general setting of a locally finite and diffuse 
A, we consider the covariance measure between the vertex point process and edge length 
process generated by Y(t). To this end, define the (random) edge-length measure €y(t) of 
Y(t) by putting for bounded Borel ACM 2 

eeMaxEdges(Y(t)) 

with £(■) standing for the usual one-dimensional length. The object of our interest is the 
measure Cov(Vy( t ), £y(t)) given by 

/ (f®g)dCov(V Y (t),£Y«t)) = Cov(E^/ (Y(t)), E Jg (Y(t))) (30) 

J(R 2 ) 2 

for bounded measurable /, g : IR 2 — > R with bounded support, where J 9 denotes the 
functional J 9 (e) = J e g(x)£(dx). Similarly as in (20) we have 

A Jg (Y) = [ J 9 (L)A(dL). (31) 

J[R 2 ] 

Thus, Aj g is constant and hence, using (3) for <fr = rjf and tp = J 9 , taking expectations 
and recalling (18), yields 

Cov(E, / (F(t)),E JS (r(t)))= / EA v f J9 (Y(s))ds+ 

Jo 

2 rCov(E A/H (F( S )),E j9 (r( S )))rf S . (32) 
Jo 

Using (3) once again, with = A'[-] and = J 9 , upon taking expectations and recalling 
(31) and (20), we get 

Cov(E A/N (y( s )),S j9 (y( s ))) = [ S EA Af[ . ]Jg (Y(u))du. 

Jo 
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Substituting into (32) leads us to 

Cov(S T;/ (F(t)),S j9 (F(t))) = / EA vfJg (Y(s))ds 

Jo 

+2 f [ EA Af[ . ]Jg (Y(u))duds. 
Jo Jo 

Applying (24), we obtain therefore 

Cov(E r? /(y(t)), Hjg(Y(t))) 
= \l #)J 9 (e)I 1 (s 2 expHA([e]));<)(((AxA)nA})(&) 

+ [ A f [e]J 9 (e) X 2 (s 2 exp(-sA([e]));t)(((A x A) n A))(de). (33) 

Consequently, putting (33) together with (30) and defining the measure (£(• He))(A) := 
e(Ane), ACM 2 , we obtain 

Theorem 3 For general locally finite and diffuse driving measures A we have 
Cov(Vy {t , w) ,£ Y{t , w) ) = f (l-A e ®£(-ne) X 1 ( S 2 exp(- S A([e]));t) 

J[R 2 ] \ Z 

+A[- n e] ® £(■ n e) X 2 (s 2 exp(-sA([e])); t)) (((A x A) n A))(de). (34) 

Observe, that the first term in the integral (34) takes into account the pairs consisting of a 
vertex constituting the terminus of a maximal edge and the maximal edge itself, whereas 
the second term corresponds to pairs consisting of a vertex lying in the relative interior of 
a maximal edge and the maximal edge. In analogy to the case of Theorem 2, vertex-edge 
pairs where the vertex is not adjacent to the edge bring no contribution to the considered 
covariance structure. 

3.4 Variance Calculation for the Stationary and Isotropic Case 

We specialize now the results obtained in the preceding Section 3.1 to the stationary and 
isotropic case, i.e. we consider stationary and isotropic random tessellations in the plane 
that are stable under iteration (STIT tessellations). Up to reparametrization, it means 
taking the driving measure Aj SO - the isometry invariant measure on the space of lines in 
the plane with length density one. Recall first that Crofton's formula [11, Thm. 5.1.1] for 
d = 2 and k — 1 implies 

A iso ([K}) = -Vi{K) = -P(K), (35) 
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where if C I 2 is a planar convex body with first intrinsic volume V\(K) and perimeter 
length P(K). For a line segment eel 2 this is just 

Kso([e}) = -£(e), (36) 

71 

where £(e) stands for the length of e. It follows from (35) that 

Z Atso[ . ] (Y(s,W)) = h e (Y(t,W)). 

Moreover, in the context of (6) and (7) we have 

2 2 

{{A iso nA iso )){dx) = -dx and A K . so[ . ] {Y(s, W)) = - Area(W). 

7T 7T 

Recall from Theorem 4 in [12] that the variance of the total edge length in W of the 
stationary and isotropic iteration stable random tessellation Y(u, W), u > simplifies in 
our particular case to 

/°° / 9 \ dr 

Iwir) (l-e~ w ) - 

where 7jy(r) = f s Are&(W D (W + ru))vi{du) is the isotropized set-covariance function 
of the window W, with v 1 standing for the uniform distribution on the unit circle «Si, see 
[13] for the definition of 'Jw(') an d Subsection 4.2 in [12] for further details. Combining 
this with (8) and (12), in view of (36) we are immediately led to the variance formula. 
An alternative method for deriving this formula directly from (16) in Theorem 1 in the 
case A = A iso is to use (35) and the important identity for the intersection measure (15), 
namely 

{{(Aiso x A lso ) n A lso )){dxy) = f X V (37) 

7r 6 \\x — y\\ 

The last equation may be established by a twofold application of the afline Blaschke- 
Petkantschin formula [11, Thm. 7.2.7] as shown in [12, Eq. (50)] in connection with (36). 
It follows 

= <P W + ^Area W 2 -2/ /gH^O ^ 
vr JwJw (l\\x-y\\) 3 ^\\x-y\\ 

= -P(W) + -Aiea(W)t 2 -2tt / i w {r)— — \ " > dr. 

71 TV JO r 

Computing now T| xp (— |ir), we arrive at 
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Corollary 1 The variance of the number of maximal edges of a stationary and isotropic 
random iteration stable tessellation Y(t, W) is given by 

Var(E,(y(t,H0)) = 



*-P(W) + \ Area(H--)« 2 + jf ' J w (r) - * + H* (l - «"*")) *. (38) 

As an example we may consider for the ball B 2 R in R 2 with radius R > 0. In this special 
case, the isotropized set covariance function takes the special form 

7 B a (r) = 2i? 2 arccos - ^v^i? 2 - r 2 , < r < 2R. 

Unfortunately, the arising integral cannot further be simplified. 



We are now interested in the variance asymptotics for a sequence W R = R ■ W of growing 
observation windows, as R — > oo. To this end, first note that asymptotically, as R — > oo, 
we have 



dr ~ ^ 2 logi? 



as long as log-B(-R) ~ logi?, where -B(-R) stands for some upper integration bound depend- 

f(R) 

ing on R. Here and later, we will write f(R) ~ whenever lim = 1. Now, the 

r^oo g{R) 

relation 

7V fl ~ Area(W^) = i? 2 Area(Vy), 

valid uniformly for the argument r = 0(R/ logi?) and 7 W — > for r = f2(i?logi?) (using 
the standard Landau notation), implies 

Var(Ei(y(*,W fl ))) ~ -tP(W R ) + -t 2 km&(W R ) + -t 2 A.m&{W R ) log R 

71 71 71 

~ -t 2 Are&(W)R 2 logR. 

71 

Summarizing, we have shown 

Corollary 2 Asymptotically, as R — > oo, we have 

Vax(Ei(Y(f,WB))) ~ -Area(W)t 2 i2 2 logi2 

7T 

and 

Var(A^(Y(i, W^))) Area(W)t 2 R 2 log R 

where N v (Y(t, W R )) ~ 2E 1 (Y(t, W R )) is the number of vertices ofY(t, W R ). 
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Figure 2: Variance comparison for the number of vertices of Poisson line (PLT), STIT and 
Poisson-Voronoi (PVT) tessellations for t — 1 and Wr = B R) the ball with radius R > 



The formulas show that the geometry of the window W is only reflected by its area in the 
variance asymptotics. Recall from [4] that for the Poisson-Voronoi tessellation PVT(t, W) 
restricted to some window W C M 2 with edge-length density t > (i.e. the intensity of 
the underlying Poisson point process equals ^-) we asymptotically have 

Var(A,,(PVT(t, W))) ~ 2t 2 R 2 Area(W). 

Moreover, for the stationary and isotropic Poisson line tessellation PLT(t, W) with length 
intensity t and restricted to W we have according to [3] 

Var(^(PLT(t, W))) ~ ^t 3 R 3 CPI(W; 2), 

where CPI(W; 2) is the 2-nd-order chord-power integral of W in the sense of [11, Chap. 
8.6]. The appearance of CPl(W;2) in the latter formula means that beside the area of 
W, also its shape plays asymptotically an important role. The variances of the different 
tessellation models are illustrated in Figure 2 for a a sequence Wr = B R of circles with 
radius R. In this case it holds CPl(B 2 ; 2) = ^ as can be concluded from Thm. 8.6.6 in 
[11] with a corrected constant. 

In particular, the formulas from the last corollary establish weak long range dependencies 
(cf. [10]) for the point process of maximal edge-midpoints and the point process of vertices, 

S1I1CG 

,. V^i(Y(t,W R ))) Vax(N v (Y(t,W R ))) 

llm 1 777T^ = llm « 777T-\ = °°- 

R-^oo Area(n-R) R^oo Area(VvR) 
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As explained at the end of Subsection 3.2, these dependencies are propagated by long 
maximal edges on which the vertices are lying and the log-term in the asymptotic vari- 
ance formula reflects the weakness of these long range dependencies, whereas in con- 
trast to the STIT model, Poisson line tessellations have strong long range dependencies, 
since Area(W / j?)~ 1 Var(A^(PLT(t, Wr))) grows polynomially (linearly) in R, while Poisson- 
Voronoi tessellations do not have long range dependencies at all. In fact, the maximal 
edges almost surely have finite length and there are no full straight lines in the tessellation 
Y{t). In contrast to this, Poisson line tessellations consist by definition only of full lines 
and the spatial dependencies are in this case much stronger due to the geometric struc- 
ture of these processes. For Poisson- Voronoi tessellations, we have independence of local 
geometries whenever the observation regions are far enough from each other, which means 
that we have 'almost' independence for the point process of vertices, see [4]. In this sense 
the STIT tessellations exhibit features intermediate between Poisson- Voronoi and Poisson 
line tessellations. 



3.5 Vertex Pair- Correlations for the Stationary and Isotropic 

Case 

Having made in Section 3.2 very general computations for the covariance measure of the 
point process of vertices of an iteration infinitely divisible random tessellation Y(t, W), 
we specialize now - as in the last Subsection 3.4 - to the stationary and isotropic set- 
up and consider a stationary and isotropic random tessellation Y(t, W) that is iteration 
stable, i.e. a stationary and isotropic random STIT tessellation with driving measure Aj SC . 
Firstly, recall the relation (37). Plugging this expression into (29), using in addition that 
Aj SO ([- He]) = |£(- fl e) and applying the substitution y = x + u yields 

Cov(W (t) )=/ / ^A^^A^o^^^^e-^l+jt) 
Jr 2 Jr 2 

+ ( a^™ <g) -e(-nOu) + -£(■ n Ou) <g> a°" j o (^y 1 z 2 ( s V- IH|s ;t) 

\ n n J 

+ 4 (-£{■ nO^)® -£{■ n Ou)) o (tf^)- 1 z 3 ( s 2 e -^ IM|s ; t) 

\7T 7T J 



Adxdu 



7T 3 \\u\ 



where we have abbreviated by O the origin and the usual length measure by I and where i? x 
stands for the diagonal shift $ x (v , w) — (v + x, w + x), v, w e 1R 2 . The covariance measure 
Cov(Vy(i)) can be reduced in the sense of [2, Sec. 8.1] and, by Proposition 8.1.1(b) there, 
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the reduced covariance measure Cov(Vy(t)) has the form 

"1 



Cov(TV(t)) 



±(5 u + 5 _ u + 28 ) X\s 2 e-^W s -t) 



+ -£(■ n (-u)u) X 2 (s 2 e-- I|u||s ;t) 

7T 

+4 {- 2 [_ f_5 v _ w £(dv)£(dw)) l^e-^-t) 1 1(1 " 

\ n JouJOu J 



7T 3 \\u\ 



where we have used the fact that J 5 V - W (A° U <8> A° u )d(v,w) = 5 U + 5_ u + 25 G and 
/ 8 v - w (A° u ®£(-C\Ou))d(v , w) = 2£(- n (— m)-u). Recall now, see again [2], that the measure 
Cov (Vy(t)) and the reduced second moment measure /C(Vy(t)) are related by 

C^v(Vy (t) )=/C(Vy W )-A 2 4 2 , 

see [2, Eq. (8.1.6)], where £^? is the Lebesgue measure in the plane and where A stands for 
the intensity of Vy(t)- Hence, taking into account that the vertex intensity A equals H 2 , 
see [9], and transforming into polar coordinates gives us 

a n2ir poo o 

/C(Vy (t) ) = — / / (5 + 6 reiv )l 1 (s 2 e-^ s -,t) + -£(-nOre i 'f)l 2 (s 2 e-^ s -t) 
71 Jo Jo 71 

+ 4 j\ r - P)<W V) X 3 (s 2 e-^ s ; i)drd V + Qt 2 ) £ 2 R . 
From the last expression we can now calculate Ripley's K-function 

often also considered in the factorial version K(R) with 

K{R) = k{R)+{^)\m), 
see [2, Eq. (8.1.12)] or [13, Chap. 4.5]. We obtain 

K(R) = 4 / (1 + l[r < i?])X 1 (s 2 e-^ rs ; f) + - min(r, R)X 2 (s 2 e-^ rs ; t) 
t Jo n 

32 ( 1 \ 2 

+ — ( r min(r, fl) - - min(r, R) 2 J Z 3 (sV- rs ; t)dr + 7ri? 2 . 
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Splitting the integral into two parts, one integral over [0, R] and another over [R, oo), yields 



K(R) = nR 2 + - A f R 2l\s 2 e^ rs ;t) + -rl 2 (s 2 e~^ s ;t) (39) 
* Jo 71 



+ ^r 2 l 3 (s 2 e-^ s ;t)dr 

TV 2 



+^ I 1 ^";*) + -Rl 2 (s 2 e^ rs ;t) 

t J R 7T 



+p [rR-^jX\s 2 e-^-t)dr. 

Using (26), (27), (28) we finally obtain by using the definition 

g (r) = J-^-K{r) = -L-^K(r) 
yv ; 2nrdr w 2vrr dr w 

of the pair-correlation function, describing the normalized vertex density in the distance r 
from a typical vertex, 

Corollary 3 The pair- correlation function of the vertex point process Vy(t) of a stationary 
and isotropic random STIT tessellation Y(t) with edge-length density t > equals 

( \ _ 2 71 71-2 Z' 1 71 71-2 ^ ~-tr fAffl 

Note in this context that in the case of a stationary and isotropic Poisson line tessellation 
with intensity t > 0, the pair-correlation function of the point process of vertices takes the 
form 

gPLT(t) {r) = 1+ 4 

ixtr 

which can easily be concluded from Slivnyak's theorem for Poisson point processes, see [13] 
(here applied to the Poisson process of lines). A comparison of the pair-correlation functions 
g(r) of the STIT tessellation and g PLT ^>{r) and that of a Poisson- Voronoi tessellation 
is shown in Figure 3. However, in contrast to Poisson line and STIT tessellation, the 
structure of the pair-correlation function of the point process of vertices of a Poisson- 
Voronoi tessellation is much more complicated. It may be expressed by a sum of integrals 
of rather involved functions, which cannot be explicitly evaluated. For details and the 
non-trivial numerical computations we refer to [4]. 

It is interesting to compare the pair-correlation formula in Corollary 3 with the information 
on the variance of the number of edges provided in Corollary 1. To this end, use the variance 
formula 

(1 \ 2 f°° 2 
V&r(N v (Y(t,W))) = 2tt i^-t 2 J J 7 w (r)[g(r) - l]rdr + -t 2 Area(W), 
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Figure 3: Pair-Correlation function of the point process of vertices of a Poisson line (PLT), 
STIT and Poisson- Voronoi tessellation (PVT) (left) and the cross-correlation function of 
the vertex point process and the length measure (right) of a Poisson line (PLT) and a STIT 
each time with edge-length density 1 



see [13, Eq. 4.5.7] and Girling's formula thereafter, and compare it with the variance 
expression in Corollary 1. Taking into account that the number of vertices is, modulo 
boundary effects, twice the number of edges, we should have agreement of the leading 
term ^ in the integral in (38) with | • 2nr (^t 2 ) times the leading term ^72 of g{r) — 1 
in (40), where \ comes from switching between edge and vertex counts, the factor 2nr 
comes from transformation into polar coordinates and the remaining factor is the squared 
intensity of the vertex point process. Comparing these expressions we readily obtain the 
required agreement of leading terms determining the prefactor in the 0(R 2 log i?)-variance 
asymptotics. The remaining lower order o(l/r 2 ) terms in the pair correlation function 
(40) do not have to and do not agree with their suitably normalized counterparts in (38), 
because the latter contains additional area order corrections and, moreover, takes into 
account the boundary effects caused by edges hitting the boundary 8Wr of Wr. 
Another aspect that can be compared concerns the radial distribution function. For a 
stationary and isotropic random point process in the plane with intensity A > and K- 
function K{r) the radial distribution function p{r) is defined by 

Writing from now on p(r) for the radial distribution function of the point process Vy(t) of 
vertices of a stationary and isotropic STIT tessellation Y(t) with edge- length density t > 
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we can use (39) to conclude 



... 8 4vr 7T 2 (1 2tt tt 2 \ 2tr 

whereas for the radial distribution function p PLT W(r) of the vertex point process of a 
stationary and isotropic Poisson line tessellation with edge-length density t > we obtain 

p PLT W(r)=2t 2 r+ 8t 

7T 

from Slivnyak's theorem. It means that asymptotically we have 

p(r) ~ 2p PLT( * ) (r), as r ->■ oo. 

Remark 1 There are different normalizations available for the reduced second-moment 
measure JC and Ripley's K-function in the existing literature as for example [2] or [13]. We 
decided here not to normalize JC by one over the squared intensity, -p- but we normalize the 
K-function by that factor in order ensure that the pair- correlation function g(r) tends to 1 
as r — > oo. This is done to keep the formulas consistent with those from previous papers 
on STIT tessellations. This convention will also be adopted in the next subsection. 

3.6 Edge- Vertex Correlations in the Stationary and Isotropic 

Case 

Our interest here is focused on the cross-covariance measure Cov(W(t), £y(t)) of a station- 
ary and isotropic random STIT tessellation Y(t) with edge- length density t > and driving 
measure Aj SD . It describes the correlations between the stationary and isotropic random 
point process of vertices and the stationary and isotropic random length measure concen- 
trated on the edges of Y(t). The study of this measure was proposed in [14] and [15] and 
we recall some general definitions now. Let $i and $ 2 be stationary and isotropic random 
measures in M 2 with respective intensities Ai > and A 2 > 0. For a Borel set 5 C I 2 we 
introduce the random measures 

K 12 (B):=e[ ^(B + x^dx), JC 21 (B):=e[ ^(B + x)$ 2 (dx). 
J[o,i] 2 J[o,i] 2 

The measure /C12 describes <3> 2 as seen from the typical point of $1 and /C 2 i describes the 
measure $1 regarded from the typical point of $ 2 in the sense of Palm distributions. In 
[15] it was shown that K\ 2 {B) = K, 2 \{—B), which in particular implies for the ball 5 2 with 
radius r > the identity /Ci 2 (£> 2 ) = /C 2 i(£> 2 ). The cross- K-function K 12 {r) = K 2i (r) may 
now be introduced as 

AiA 2 
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and the cross-correlation function <?i2(r) of the random measures $1 and $2 is defined by 

1 dK 12 (r) 



0i2(r) : = 



27rr dr 



(41) 



compare with [14] an [15]. Informally, we could say that XiK 2 i(r) or A2-^i2(^) is the 
expectation of $i(-B 2 + x) or ^ 2 (B 2 + x ) a ^ the typical point x of $ 2 or $1, respectively. 
The theory is now applied to our setting and we take for $1 the point process Vy(t) and 
for $2 the random measure £y(t) (recall the definitions from Section 3.3). We have in this 
special situation Ai = -t 2 and A 2 = t, see [9]. We can now use Theorem 3 together with 
(37) and (35) to obtain, under the substitution y = x + u, 



Cov(Vy( t ),£w t) ) 



2 



nOu)T(s 2 e--^ s ;t) 



7T 



_ x \.dudx 

7T 3 I \u\ I 



Consequently, we end up with 



12 



+ 



- / f_5 v . w £(dv)£(dw)l 2 (s 2 e'^ s ;t) 

71 Jou JOu 



<Ou JOu 

whereby, upon passing to polar coordinates, we have 

+ 



Adudx 2 3 

on || + —t %2, 
7T \U\\ 7T 



K 



12 



jo 



(■nOre^)l\s 2 e~* rs ;t) 
~k (/ ^ ~ P^ 5 P elvd p) ^Ve - *™! t)d<pdr + ^t 3 4 



Recalling the definition of if 12 and using again that in our set-up Ai = M 2 and A 2 = t, we 
obtain 

4 Z" 00 2 
K 12 (i?) = 7rR 2 + — j mm(r,R)I 1 (s 2 e~ rs ;t) 

+^ min(r, i?) - ^ min(r, i?) 2 ^ Z 2 (sV- rs ; i)dr 
= ttR 2 + ^- f R ' rl\s 2 e-$ rs ;t) + -r 2 l 2 (s 2 e-^ rs ;t)dr 

Tit 6 Jq 7T 

4 Z" 00 



; i) + - [rR - — X 2 (s 2 e— rs ; t)dr. 



Using now (41) together with (26), (27) and (28) we arrive at 
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Corollary 4 The cross- correlation function of the vertex point process Vy(t) and the ran- 
dom length measure Sy(t) of a stationary and isotropic random STIT tessellation Y{t) with 
edge-length density t > equals 



#12(0 1 + fV ^3^,3 y 2t 2 r 2 4t 3 r 3 J e 

In contrast to this formula, the same cross-correlation function g 12 (r) for a stationary 
and isotropic Poisson line tessellation with edge-length density t > is given by 

which can easily be obtained from Slivnyak's theorem for which we refer to [13]. A com- 
parison of both functions is shown in Figure 3. We would like to point out that the 
corresponding cross-correlation function is unknown until know for the Poisson- Voronoi 
model. 



4 Central Limit Theory 

In this section we will study the functional central limit problem for the total edge count and 
edge length processes induced by a STIT tessellation in growing windows Wr = RW, R — > 
00, with W standing for some compact convex set of non-empty interior, to remain fixed 
throughout the section. We assume that the measure A is translation invariant, i.e. we are 
in the STIT regime. 

Define the rescaled total edge length process 
as well as the rescaled total edge count process 

^(y^ + i/iog^Wfl)), * e [0,1]. 



Ry/hiR 

The main result of this section is 

Theorem 4 The processes (Cf ,w , Cf : ' w/ ) sg [ 0j i] converges jointly in law, as R — > 00, on the 
space T>([0, 1]; M 2 ) ofM, 2 -valued cadlag functions on [0, 1] endowed with the usual Skorokhod 
Ji-topology, [5, Chap. VI. 1] or [1, Chap. 3, Sec. 14], to the process t 1— > where £ 

is a normal random variable with variance V(A[-],W) which is given by (63) in [12] for 
4> = A[-] there. 

Rather than giving a general formula for V((f>, W), we refer the reader to Proposition 1 and 
(72) of [12], where the general case is considered and V((f>, W) is expressed as a weighted 
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mean width of an associated zonoid and its polar body. Here, we only mention the fact 
that for the particular isotropic case we simply have 

V(A iao [-],W) = -Aie&(W), 

71 

see the discussion following Proposition 1 in [12]. 

The phenomenon observed in Theorem 4 above deserves a short discussion. Namely, al- 
though both E A r.i(F(t, Wr)) and Ei(Y(t, Wr)) exhibit fluctuations of the order RlogR, 
the mechanisms in which these fluctuations arise are of a rather different nature: 

• As shown in [12, Thm. 6], the leading-order deviations of E A [.] arise very early in 
the course of the MNW-construction, in its initial stages usually referred to as the 
big bang phase. Here, this is the time period [0, 1/ logi?]. During the later stages of 
the construction, i.e. the time interval (l/log/2, 1], the variance increase is of lower 
order and any newly arising fluctuations are negligible compared to those originating 
from the big bang. In the asymptotic picture this means the initial fluctuation remains 
frozen throughout the rest of the dynamics, whence the constant limit for C^' w (note 
at this point that the Brownian limit for the length process obtained in Theorem 6 
in [12] referred to a different time flow). 

• In contrast, the deviations of Si arise and cumulate constantly in time t with rate 
proportional to t times the initial big bang fluctuation of S A [.j. Thus, as opposed to 
that of Sa[.j, the variance of Si exhibits a non- vanishing quadratic dependency on t 
even in large R asymptotics. 

Thus, in large R asymptotics, we have the following intuitive picture: denoting by S the 
initial big bang fluctuation of E A [.j we can effectively use the following first-order approxi- 
mations: S A[ .](T(t, W R )) « S and Ei (y(t, W R )) « tE valid for t E (1/log.R, 1]. 

Proof of Theorem 4 Consider the auxiliary process 

Ei(y(£, W)) :— Ei(y(£, W)) - [ A 1 (Y(s, W))ds, 

Jo 

which is in view of (5) the same as 

Si(Y(i, W)) — [ EAuiY&Wfids 
Jo 

and which is a centered S r martingale by (2). Squaring and taking expectations we get 

E(Ei(F(t,W0)) 2 = Vax(E 1 (Y(t,W)))-2 f E[Ei(y(*, W))E A[ .](y( S , W))]ds 

Jo 

+ e[ [ E A[ . ] (r( S ,H/))E A[ . ] (F( M ,H/))^ S . (42) 
Jo Jo 



25 



Using that for s < t 

E(Z A[ . ] (Y(t,W))\%) = Z A[ . ] (Y(s,W)) 

and 

E$ 1 (Y(t,W))\S a )-Z 1 (y(s,W)) = 

^(^A[-](Y(u, W))\$t s )du — (t — s)T, A [.](Y(s, W)), 

as follows by the martingale property of T*\[.](Y(t,W)) and Tii(Y(t,W)), we can rewrite 
(42) by taking first conditional expectations as 

E(S!(y(f,wO)) 2 

= Vax(Ei(y(*,W)))-2 f (t - s)EE\u(Y(s,W))ds 

Jo 

- 2 [ EE 1 (Y(s,W))£ A [.](Y(s,W))ds + 2 [ [ ET, 2 Al] (Y(u, W))duds 
Jo Jo Jo 

and hence, with the second and fourth term in the right hand side canceling out, 

E(t 1 (Y(t,w))) 2 = 

Var(E 1 (y(t, W))) - 2 / Cov(E 1 (y(s, W))Y. m {Y{s, W)))ds. (43) 

Jo 

The relation (43) combined with (10) and (8) readily yields 

E(E 1 (y(t, W))f = tA([W}) + ^«A n A))(W). (44) 
To proceed, define now the auxiliary process 

C? W := —±==£ 1 (Y(t + l/\ogR,W R ))=C?> W - f Cf' w ds. (45) 
RVlogR J-i/iogR 

Using that A([W R }) = 0(R) and ((A n A))(W fl ) = 0(R 2 ) we conclude from (44) that 

lim E(C?' W ) 2 = 0. 

R^O 

Thus, Doob's L 2 -maximal inequality [6, Thm 3.8(iv)] implies 

lim E sup (C*' w ) 2 < lim 4E(Cf ,w ) 2 = 0. (46) 

We are now in a position to apply Theorem 6 in [12] (taking into account that the notation 
£R,w usec } in [12] corresponds to a time change of L R)W as defined here) to conclude that, 
as R — > oo, 

{Cf )W )te[o,i\ =>• (0*6[o,i]> 
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that is to say the process (C t ' )te[o,i] converges in law in the Skorokhod space V([0, 1];M) 
to the constant process i 4 where, recall, £ is a centered normal random variable 
with variance V(A[-], W). Using (46) and recalling the definition (45) we see now that the 
processes (Cf' W ,C t ,W )te[o,i] converge jointly in law in T>([0, 1];R 2 ) to the process 



Remark 2 For reasons discussed in detail in Remark 6 in [12], we expect the rate of the 
convergence in Theorem 4 to be rather slow. 
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